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Abstract 

We present the general solutions for the classical and quantum dynam- 
ics of the anharmonic oscillator coupled to a purely diffusive environment. 
In both cases, these solutions are obtained by the application of the Baker- 
Campbell-Hausdorff (BCH) formulas to expand the evolution operator in 
an ordered product of exponentials. Moreover, we obtain an expression 
for the Wigner function in the quantum version of the problem. We ob- 
serve that the role played by diffusion is to reduce or to attenuate the 
the characteristic quantum effects yielded by the nonlinearity, as the ap- 
pearance of coherent superpositions of quantum states (Schrodinger cat 
states) and revivals. 



PACS: 03.65.Yz: Decoherence; open systems; quantum statistical meth- 
ods, 02.20.Sv: Lie algebras of Lie groups 



1 Introduction 

In the last dècades, the investigation about the transition from quantum to 
classical dynamics has progressed enormously. This was induced in part by the 
development of the experimental techniques, especially in quantum òptics pQ, 
and in part by the possibility of appearance of technology in quantum infor- 
mation processing [2], (3]. Such research aims to understand how the typically 
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quantum efïects disappear in the dynamics of macroscopic systems. According 
to a popular theoretical model, onc believes triat the cmcrgcncc of the classical 
world from quantum mechanics is a consequence of the unavoidable coupling 
between the macroscopic system and its environment. For example, in accor- 
dance with this proposal [3] , environmental coupling is responsible by the ràpid 
evolution of coherent quantum superpositions of macroscopically distinguishable 
states (Schròdinger cat states) into statiscal mixtures, a phenomenon known as 
decoherence. 

Despite these theoretical advances, some aspeets of the quantum to classical 
transition remain subtle and controversial, especially in the case of classically 
nonlinear or chaotic systems. For some authors, the departure of the quantum 
mean vàlues of observables from the corresponding classical ones (correspon- 
dence breakdown) in chaotic systems oceurs in a very small time scale (see, 
for example, [5]) and they sustain the idea that the coupling with purely diffu- 
sive environment can reduce these discrepàncies and increase the break time by 
decoherence. On the other hand, for other authors (see [ïï]) decoherence is not 
necessary to explain the classical behavior of macroscopic systems (including the 
chaotic ones) since the observed discrepàncies between quantum and classical 
mean vàlues of observables are negligible for any current realistic measurement. 

In order to shed some light on this debaté, we turn our attention to the 
one-dimensional nonlinear or anharmonic oscillator (AHO) coupled to a purely 
diffusive reservoir. Three are the reasons for the choice of this model. Firstly, in 
the limit of vanishing environmental coupling, its quantum dynamical evolution 
exhibits several effects without analogous in its classical counterpart, such as 
revivals and appearance of coherent superpositions of states [7]. Secondly, by 
virtue of its relative simplicity, it is possible to obtain the exact solutions of the 
equations of motion for quantum [HI El HHl HH E^l H31 E3] and classical versions of 
the model even in the presence of the reservoir. Last, but not least, the recent 
technical advances in trapping and controlling cold àtoms suggest the Bose- 
Einstein condensates (BEC) as potential candidates to implement experimental 
tests of this model. In fact, in the single-mode approximation, a BEC trapped 
in a optical lattice is suitably described by the quantum nonlinear oscillator 
[15] . Hence, dissipative AHO became a largely studied model in the literature. 
In Ref . [TD] , Daniel and Milburn obtained the exact evolution of the Q func- 
tion associatcd to an initial coherent state of an AHO subject to attenuation 
or amplification. The authors showed that the effects yielded by nonlinearity, 
such as revivals and squeezing, gradually vanish if a non-unitary mechanism is 
taken into account. These results were extended by Pefinovà and Luks [TTJ to 
an arbitrary initial state. In Ref. [13], Kheruntsyan obtained the steady state of 
Wigner function of a single driven damped cavity mode in the presence of a Kerr 
medium. The author studied a model of reservoir that included two-photon ab- 
sortion, besides the usual one-photon absorption. Closely related to the present 
contribution, Chaturvedi and Srinivasan [12j found an exact solution of a class 
of master equations governing the dynamics of a chain of coupled dissipative 
AHO. The authors used thermofield dynamics notation in order to map a mas- 
ter equation into a Schròdinger equation with a non-hermitean hamiltonian. 
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The form of the general solution of this equation was expressed as an ordered 
product of exponentials of operators acting on an arbitrary initial condition. 

In this contribution, we obtain the exact time evolutions of the density op- 
erator (section 2) and of the classical distribution function (section 3) in the 
quantum and classical versions of the nonlinear oscillator coupled to a purely 
diffusive environment. Both results are obtained by the application of the Lie 
algebraic techniques, in particular, the Baker-Campbell-Hausdorff (BCH) for- 
mulas used to expand a Lie exponcntial in an ordered product of exponentials 
[T2I ÏÏM [TTl ÍT51 [TÜ] . In the quantum case, the result allows us to find the corre- 
sponding Wigner function in terms of an expansion in associated Laguerre poly- 
nomials. It is important to mention that the Wigner function should smoothly 
approach the corresponding classical distribution in the appropriate limit [5|. 
Thus, we take the classical limit of the partial differential equation (PDE) that 
governs the time evolution of the Wigner function and we obtain a classical 
Fokker-Planck (FP) equation. The algebraic structure of the master equation 
in the quantum case is preserved by the corresponding FP equation in the clas- 
sical case, so that we can extend the methods applied to obtain the solution of 
the first in the finding of the solution of the last. 

In Rcf. 20J, Oliveira and co-workers investigated the diffusive AHO and 
showed that the break times (i.e., the characteristic times of departure of the 
quantum and classical dynamics) depend strongly on observable and initial con- 
dition. A more "fair" comparison between the two dynamics shall be given by 
the evaluation of the distance between the corresponding distributions in the 
phase space [2T]. Hence, the exact solutions of the equations of motion for the 
Wigner function and classical distribution function will allow to define analyt- 
ically the break time for the diffusive AHO and its dependence in terms of the 
nonlinearity strength and diffusion constant. 

We conclude this work by comparing the quantum and classical evolutions of 
the Wigner and the classical distribution functions for an initial coherent state 
in the AHO with and without diffusion (section 4). The results are preliminar 
and deserve a more careful analysis, but they suggest that, as expected, in 
quantum case, inclusion of diffusion reduces the phase space interferences and 
therefore prevent the appearance of quantum coherent superpositions. The area 
of the regions where the Wigner function should be negative is reduced too. As 
time goes by, the Wigner distribution gradually takes the form of an annular 
volumc around the origin of the phase space. For later times, this volume 
becomes more "fat" and "flat" . In the classical case, the fine-structured whorl 
yielded by the distribution in the diffusionless regime is destroyed. As it happens 
in the quantum case, diffusion turns the classical distribution more "fat" and 
"flat" around the origin. These results suggest that the Wigner function of the 
quantum diffusive AHO converges gradually to the distribution function of the 
corresponding classical version of the model in non-unitary evolution. 
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2 Quantum mechanical diffusive anharmonic os- 
cillator 



Let us consider the AHO coupled to a thermal bath of oscillators in equilibrium 
at temperature T. Assuming that the nonlinearity strength is small and the 
coupling to the reservoir degrees of freedom is weak, we obtain the following 
master equation in the interaction picture: 

m = (i) 



-w 



(à f à) 2 ,p(i) + k(n + l) [2àp (t) à f - iïàp (t) - p (t) a f fi] 



+kn Ï2à' p (i) à — aaJ p (t) — p (t) àà'] , 

where k and g are the damping and nonlinearity constants, respectively, and 
fi is the average number of thermal photons in the mode lo of reservoir (uj is 
the natural frequency of the oscillator) . The density operator p (t) represents 
the state of the system at time í; à and à> are the annihilation and creation 
operators, respectively. We are interested in the so called diffusive limit of the 
above equation. This limit is obtained by taking the damping constant going 
to zero, k — > 0, and the number of thermal photons going to infmite, n — > oo, 
keeping the product k — kn finite. Thus, the master equation (fï} becomes 

p(t) = C^pit) (2) 



-i<) 



(a f a) ,p(t) 



+2k [à/5 (í) a/ + a 1 p (t) à — à^àp (t) — p (t) a'à — p (í)l 
= -ig (M 2 - V 2 ) p (t) + 2k [J + K - {M + V + 1)] p (t) . 

Using the notation given in Ref. [55], the super-operators (sup-op) that 
appear in the last line of above equation are defined by 

M = (a^à), = àjài , (3a) 

V = (à f à) r = à r àl , (3b) 

J = àià\ = à\ài , (3c) 

1Z = à r à\ = à\à r . (3d) 

Here, ài, fií, fi r , fij represent the left and right actions of the creation and 
annihilation operators on a genèric operator O: 



ai 



Ó = àÓ, à\Ò = à^Ò , a r Ò = Òà , atÓ = Ófit . (4) 



At this point, it is convenient to introduce the following nomcnclature. We 
call B (H) the set of operators that act on the oscillator space state TL. The 
elements of B (H) can be assigned to vectors of an extended Hilbert space con- 
structed by direct product between the original space state and its dual 7í*, viz. 
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7í<S>'H*. This extended Hilbert space is frequcntly called Hilbcrt-Schmidt space 
or Liouville space [23] , 

The formal solution of Eq. © is given by 

p(t)=e c °° t p(Q) , (5) 

where /5(0) represents the initial state of the AHO. The evolution of a genèric 
initial state 

P(°) =^2Pm,n\m) (n\ (6) 

can be evaluated by expanding the exponential e Coat in an ordered product of 
cxponcntials. Usually, this task is achieved by the systematical application of 
the Lie algebraic methods, in particular, using the Baker-Campbell-Hausdorff 
(BCH) formulas PIIQIIIIIQIICI]. For the dissipative AHO, this expansion was 
obtained in Ref. [12] . and we reproduce in detail the procedure here. To carry 
out this expansion, we begin by evaluating the commutation relations between 
the sup-op defined in (J3j) . They are listed in Table [U 





M-V 


\ (M + V - 


- 1) n 


J 


M -V 














\(M + V + 1) 








K 


-J 


K 





—11 





-(M+V + l) 


J 





J 


M + V + l 






Table 1: Commutation relations between the sup-op defined in Eq. Q. The 
i,j entry in the table is the result of the commutation between the sup-op of 
the í-th row and the sup-op of the j-th column. 



The sup-op M —V, \ (Ai + V + 1), J, 1Z form a four-dimensional Lie àl- 
gebra, which we will denominate A4. We easily recognize a subalgebra su (1,1) 
contained in A4, formed by the set \ (M. + V + 1), J, 1Z. Let us rewritten the 
Liouvillian as 

Coo = ig (M - V) - A (M + V + 1) + 2k ( J + TZ) , (7) 

where we define the sup-op 

A = ig (M -V) + 2k, (8) 

that is to be formally considered a c-number, since the sup-op M — V commu- 
tates with the rest. The eigenvectors of A are the same of M — V. It is easy 
to verify that these eigenvectors belong to the set {|m) (n\}, with eigenvalues 
ig (m — n) + 2n. 

The formal solution ([5]) can be written as 

p (t) = exp [-Aí (M + V + 1) + 2nt (J + K)} exp [igt (M - V)} p (0) . (9) 
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Applying the well-known BCH formulas to expand su (1,1) Lie exponentials 
[Ï9l [24l [25] , we rewritten the above expression as 



p(t) = exp [r (t) TZ] cxp {ln [r (í)] [M +V + 1)} exp [r (í) J] 

xexp[igt{M -V)]p(0) , (10) 

where T (t) , and Tq (t) are time-dependent sup-op given by 

r ° (í) = A cosh (Ai) + Asinh (Ai) ' 
. . 2k sinh (Aí) 

= A cosh (Aí) + Asinh (Ai) ' ( ' 

Here, we define 



A = VA 2 -4k 2 . 
Let us consider the initial state (J6j) . Its time evolution is 

p(t) = exp[r(í)^]exp{ln[r (í)](M+'P + l)}exp[r(í) l 7] 
x pmn exp [igt (to - n)] |m) (n| . 

The action of exp [r (í) ^7] on |to) (n| is obtained in the following way: 

J2^p[^(t)j}\m)(n\ = j2T,^r a3 \ m ^ n \^y ( 12 ) 



^ / r J (í) 

where (m, n) = min (m, n), and 7„ (í) is a c-number function obtained by direct 
substitution of the sup-op A by ign + 2n in the expression (jllb[) . Let us redefine 
the indexes m — > m — j , n — > n — j . Thus, 

/3(í) = exp[r(í)^]exp{ln[r (í)](X+7' + l)} (13) 



v— / (to + j) \ (n + j): 
x 2^ 2^ Pm+j,n+j y — j-j exp [ï 5 í (to - n)J 

j m,n 

x -| M («I ■ 

The action of exp {ln [r (i)] (M. + V + 1)} on |to) (n| is evaluated in ana- 
logue way: 



exp {ln [To (t)] (M + P + 1)} |m) (n| = T (i) £ {ln [F ° (t)]}3 (X + V) J |m> (n\ . 
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\m) (n\ is an eigenstate of Ai + V with eigenvalue (m + n). Hence, 
exp{ln[T (t)](M + V + l)}\m){n\ = r (t) ^ ^ ( m + n ) j \m) (n\ 

= rZ +n+1 (t)\m)(n\=C±: +1 (t)\m)(n\ . 

Here, £ n (t) is a c-number function obtained by the direct substitution of the 
sup-op A by ign + 2k in the expression (| 1 1 a,|) . Substituting this result on Eq. 
(fT3|) . we have 

p (t) = exp [r (t) 72.] 



2_, ftM-i,n+j V -jZl ex P 1*5* ( TO ~ n )J ( 14 ) 



. . . (m + j)l(n + j)\ 

x / / Pm+i.n+i\l , , 

mini 

xC +„+i (Í) 7^!)| TO) (n |. 



The action of exp [r (í) 72,] on |m) (n| is evaluated in the following way: 
exp[r(í)72-]|m)(n| = ^^-jf- (a*) 1 \m) {n\ a l 



7m-n (*) / (m + Q!(n+7j! 

> 7i V t—x \m + l)(n + l\ 

l\ V m!n! 1 ' x 1 



Substituting this result in Eq. (|T4[) , we finally have the general solution of Eq. 
([2]) with the initial condition ([6]): 



'(*) = L·L·L·^>+i 

l j m,n 

X7™ + Í„ (*) Cí^ +1 (í) exp [iflí (m - n)] |m + l) (n + l\ . (15) 

2.1 The Wigner function of the diffusive anharmonic os- 
cillator 

We can obtain a representation of operators that belong to B (H) as functions 
of the set ^"(0), i.e., functions on the phase space Í7. This representation can 
be seen as an invertible mapping between B (H) and the set T (Q) . One of these 
representations is given by the Weyl- Wigner transform [SS] [S7J 113 > defined on 
a genèric operator O as 

O (q,p) = ^ é up ' n (q - u/2\ Ò\q + u/2) du = W (d) (q,p) , (16) 

where q and p are phase space coordinates position and momentum. The Weyl- 
Wigner transform of the operator density, defined as 

W (q,p) = ^ J°° é^ h (q - u/2\ p \q + u/2) du = W (j^ (q,p) , (17) 



7 



yields a quasi-probability distribution function - the Wigner function, W (q,p). 
We can merge the phase space coordinates q and p in an unique complex variable 
q as follows 

Mlü . / ï _Q + iP 
~2h + lP \l 2Muoh = y/2 ' 
Here, M is the mass of the oscillator, and Q and P are adimensional real 
variables. Adopting this representation, the Wigner function is defined by 



W{a) = \ tr \pD (a) é^ à D^ (or) 



(18) 



where D (a) = exp faà' — a*a) is the unitary displacement operator [30j . and 
tr (■) stands for trace. 

The application of the Weyl- Wigner transform in expression (fTÏÏl) produces 



W W = Prn+j.n+j 



yj{m + j)\ (n + j)\ (m + ï)\ (n + l)\ 



m!n!Z!j! 

X7„i n (*) Cí; i+1 (*) exp [iffí (m - n)] U m+l<n+l) (19) 

Hence, the Wigner function W (t) of the nonlincar oscillator is expressed in 
terms of the functions H m ,n, that are obtained by the Weyl- Wigner transform 
of the eigenfunctions {\m) (n\} m n=0 1 of the sup-op (M — V), i.e. 



n, 



(_-t\m 

(a) = (2nhy 1 W (|m) (n\) (a) = i— (n| D (2a) |m) 

7TAÏ 



The matrix elements of the operator D (2a) are given by (see appendix B of 
Ref. [30]) 



(n\D(2a) \m) 



^e" 2 !"! 2 (2a)"- m L"r m (4 |a| 2 ) 
~e" 2 H 2 (-2a*) n ~ m L™~ n (4|a| 2 ) , 



n > m 



m > n 



where L 7 ^ m (x) is the associated Laguerre polynomial. Hence, we have 



n„ 
n„ 



.(«) 



.(«) 



(-1)' 

■nh 

(-i) n 



-2|c 



.p"2|a| 



(2a) 



«— m rn— m 



(4|a| 2 ); 
(2a*) m ""L^-" (4|a| 2 ) ; 



71 > TO, 



to > n. 



(20) 



(21a) 



(21b) 



7r/ï V TO! 

We are interested in the time evolution of an initial coherent state p (0) = 
l a o) ( a o\- I n this case, the matrix elements of the density operator are p m „ = 
Oi™ (oíq)™ e^l a °l" (TO,!n!)" 5 . Substituting this into (fTÏÏl) , an d reducing the sum in 
j we have 



W(a,t) 



l«o| 2 



m!n!/ 

ï m,n 



(22) 



x exp 



|cko| j m -n (t) + igt (m - n) U m +i, n +i (a) 
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3 The classical limit of the diffusive anharmonic 
oscillator 



3.1 The equation of motion for the Wigner function and 
its classical limit 

Taking the Weyl- Wigner transform in both sides of master equation @ we 
obtain a partial differential equation for W (a). The following correspondence 
formulas EHl are useful in the execution of this task: 



ap 



à'p 



a + -d, 



1 



W {a) , pà la - -d a * ) W (a) 
l -8^j W (a) , fití -> (c 



(23) 



+ -d a )W (a) 



where <9 Q 



' and 9 Q - 



tt. The time evolution of the Wigner function for 



the diffusive AHO is governed by the PDE 



W(a,í) 



-'.9 



^2 \a\ 2 — íj (a*d a * — ad a ) — ^ (a*d a * — ad a ) d a d a 



-2Kd a d a »}W(a,t) 



(24) 



where d t = Jj. 

The classical limit of the AHO can be obtained by taking the limit H/J — > 
in Eq. (f2"4"]l . where J is a characteristic classical action. Assuming that the 
initial state is a coherent state with amplitude ao, J ~ fi |ao| • In this case, the 
classical limit corresponds to take |ao| — * oo. In order to do this, let us define a 
new phase space variable = a/ \oíq\. In terms of this variable, the PDE 
becomes 



d t W((3,t) = 



-w 



(21 



«0 



1 09* dp-pdp) 



1 



4|a | 



+2- 



K 



■dpdp*}W(/3,t) 



(25) 



\ao\ 

Here, dp = ^ and dp* = -j^p- However, the constants g e k are dehned 
in such way that this limit does not make sense, since the nonlinear term is 
proportional to |ao| 2 and the stochastic sector of this equation vanishes. In 
order to circumvent this problem, we redefine them: 



g' =g\a \ x 



and 



k' = k/ \ao\" 



(26) 



Proceeding in this way, in the classical limit, the above PDE becomes a Fokker- 
Planck (FP) equation for the classical probability distribution w (/3,í): 

dtw(fi,t) = {2ig'\l3\ 2 (pdp-f3*dp*) + 2 K 'd fj d fí *}w(p,t) , 
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Resorting to the definition of the constants g' e n' in Eq. l|2"ü)) , we have 



dtw (a, t) — ^2íg \a\ 2 (ad a — a*d a *) + 2nd a d a * j w (a, t) . (27) 

Note that the above equation contains only partial derivatives in a and a* of 
order one or two. The partial derivatives of superior order, responsible for the 
nonlocal character of Eq. ([24]) . vanish in this limit. 
At this point, we introduce the Poisson brackets 

[f,g] P = ih- 1 [(d a f) (d a ,g) - (d a g) (d a .f)] , 

and the Eq. (I27| can be rewritten as 



dtw(a,t) — —hg \a\ ,w(a,t) — 2h k [a, [a* , w (a, í)] 



ipip 



(28) 



3.2 The algebraic structure of the time evolution equation 
for the classical distribution function 

Our objective is to find the solution of the Cauchy problem described by Eq. 
(|2~T|) and the initial condition w(a,0). It is interesting to note that, whereas 
the Weyl-Wigner transform maps operators in B(7í) into functions in ^"(íí), 
sup-op are mapped in differential operators acting on The main benefit 

of this mapping is the preservation of the commutation relations. The gener- 
ality of the Lie algebraic techniques allows one to extend the results obtained 
for a problem involving a particular realization of a determined Lie àlgebra to 
another realization of the same àlgebra. However, it is necessary to remember 
that in the classical limit, this correspondence can not be complete. Compare, 
e.g., the time evolution equation for the the Wigner function (f!H| and for the 
classical distribution function (|27p . In the classical version, the nonlinear hamil- 
tonian term does not present derivatives in the phase space coordinates of order 
superior to two. 

Let us consider, for example, the sup-op M, that acts on a genèric operator 
O as follows 

MÒ = iïàÒ . 

Taking the Weyl-Wigner transform in both sides of above equation, we obtain 



W [MO] = W 



2 1 1 

\a\ - - (1 + ad a - a*d a .) - -d a d a , 



O(a). 

Proceeding in this ways, we obtain the following relations between the sup-op 
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and the differential operators: 



M - 


♦ \a 


2 - 


1 

2 


(1 


+ a<9 a 


- a*d a *) - 


^d a da- 


V - 


♦ \a 


2 - 


1 

2 


(1 


— a<9 a 


+ a*d a -)- 




J - 


♦ |a 


2 + 


1 

2 


(1 


+ ad a 


+ a*d a *) + 


^d a da- 


K - 


» |a 


2 - 


1 

2 


(1 


+ ad a 


+ a*d a *) + 


^d a da' 



(29) 



At this point, it is useful to introduce the following differential operators: 

Y"o =a*d a ' — ad a , (30a) 
Y z = \a\ 2 - ~ da da * , (30b) 

Y+ = \a\ 2 - i (1 + ad a + a*d a .) + ^d a d a - , (30c) 

Y- =\a\ 2 + ^{l + ada+a*da*) + \d a da* , (30d) 

In terms of these operators, the classical time evolution equation for the diffusivc 
AHO becomes 



d t w (a, í) 



~^(Y + + Y^+2Y Z )Y 

+2k (Y + + Y_ - 2Y Z )} w (a, t) (31) 
LecW (a,t) , 



The formal solution of the Cauchy problem is given by the application of a Lie 
exponential on the initial condition w (a, 0), 

w (a, t) = exp (L x t) w (a, 0) . (32) 

We can express the Lie exponential exp (Looí) as a product of exponentials 
of which the action on functions in T (fi) is known. For this, we use the BCH 
expansion formulas, in analogous way to the quantum version of the problem. 
The first step is to determine the Lie àlgebra generated by the commutation 
between the operators that appear in Eq. (|3l·|) . 

The operators defined in Eq. (f3"P)) obey the commutation relations presented 
in Table [2j Comparing Tables [2] and [U we easily note that the differential op- 
erators in Eq. (|30p yield another representation of the four-dimensional àlgebra 
A4. Because of this, we can identify a Lie su (1,1) subalgebra defined by the 
operators {Y z , Y + , y_}. 

In the master equation @ , the unitary nonlinear term introduces the sup-op 
M 2 - V 2 = (M — V) (M + V). The inclusion of this sup-op to the set consid- 
ered in Table [T] yields an infinite Lie àlgebra when we evaluate its commutation 
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Y 


Y z 


Y+ 


Y- 


Y 














Y z 








y+ 


-y_ 


Y+ 





-Y+ 





-2Y Z 


Y_ 





Y_ 


2Y Z 






Table 2: Commutation relations between the operators defined in Eq. (|5D|) . The 
i,j table entry is the result of the commutation of the operator in the i-th row 
with the operator in the j-th column. 



relations with another sup-op. The trick used to find the solution of Eq. ((2]) 
consists formally in considering the sup-op (Ai —V) as a c-number, since it 
commutates with the rest and the functions obtained in the expansion of the 
correpondent Lie exponential arc, in fact, functions of this operator. We can 
determino the quantum statc in timc í, p(t), by expanding the initial state po 
in terms of the eigenfunctions of (Ai —V). 

In the same way, the nonlinear hamiltonian term in Fokker-Planck equation 
(|27[) introduces produets of differential operators and they also yield an infmite 
Lie àlgebra. However, these produets are in the form Y z Yq, Y±Yo. Since lo 
commutates with the rest of the elements defined in Eq. (|30|) , we can employ 
an analogous trick to that employed in the solution of (|27|) . i.e., we can formally 
consider Y a c-number and evaluate the functions in the correspondent Lie series 
as functions of this operator. However, the solution w (a,t) determined in this 
way will be "useful" if the action of these differential operators on elements 
in T (O) is known. Instead and analogously to the procedure adopted in the 
quantum version of the problem, we prefer to find the eigenfunctions of Yq and 
to express the initial state in terms of them. If we know how each operator in 
Eq. (|5D)) acts on these eigenfunctions, the classical distribution w (a, t) can be 
written as an expansion in terms of them with time dependent coeficients. 



3.3 Eigenfunctions of Y 

We can easily determine the eigenfunctions of Yç, since they are related with the 
eigenfunctions of the sup-op (Ai — V) by the Weyl-Wigner transform. Rcmcm- 
bering, the eigenfunctions of this sup-op are {\m) (n\\ m n=0 1 , with eigenvalues 
m — n. By Eq. (|2"U|) . the Weyl-Wigner transform of the sup-op (Ai — V) yields 

= Y W (ò) , wherc Ò is a 



the differential operator Y , i.e., W (Ai -V)0 
genèric operator acting on TL. Making O = \m) (n\, we have 

W [(Ai - V) |m) (n\] = Y a W (\m) (n|) = (m - n)W (\m) (n|) . 

Therefore, the eigenfunctions of Yq are the Weyl-Wigner transform of |m) (n\, 
namely II mj „, with eigenvalues (m — n). 

Since the rest of the differential operators defined in Eq. (|2"5|) are related 
with the sup-op defined in Eq. ([3]), and the action of the last on elements of set 
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{\m) (n|} m n=0 1 is known, we directly obtain the corresponding action of the 
first on the functions Tl m ,n- In fact, we have 

5on m ,„ = (m-n) n mi „ , 

YJ\ m , n = i (m + n + 1) n mj „ , (33) 

í+n m , ri = V ( m + !)(« + i)n m+ i : „ + i . 

In the problems that we are interested, we necd to compare the time evo- 
lutions of the Wiener function and of the classical distribution associated to a 
given initial statq^j p . For our purposes, it is interesting to express such function 
as an expansion in eigenstates of Yq. Using the properties of the Weyl-Wigner 
transform W we have 

W = W (0) = W (M = P™,nW (M H) ■ 

Wg = } j Prn,nïí m ,n- (34) 



Hence, 



3.4 The time evolution of the classical distribution func- 
tion 

Consider the initial state po = J2 m n Pm,n \ m ) ( n \- K the corresponding Wigner 
function Wq, given by Eq. (I34|) , represents a vàlid classical distribution function, 
we can make wq = w (0) = Wq. The solution of Eq. (|27|) for this initial condition 
is 

w (t) = exp (Looí) w 

= exp [G (t) Y+] exp [2Y Z ln G z (i)] exp [G (t) Y_] , (35) 

where 



n í4-\ 4v%KÏo (og \ 

G z (t) = — -, -^=r -^^—r- , (36a) 

A^/ignYQ cosh (2ívigK>o) + (4/c + igYo) sinh (2t^/igKYo) 

G(t)= - sinh (2tV^yo) 

4Vi.gK^o cosh (2ty/igKY ) + (4/c + z#y ) sinh (2ty/igKY ) 

Note that G (i) and G 2 (í) are functions of operator Y . 

The action of exp [G (í) Y"_] on the initial condition wq yields 



exp [G (t) Y-]w = }^}^ Pr+j,s+j — Tj y ^ n r 



1 A classical distribution associated to a state p is a probability function that the marginal 
distributions coincide with the corresponding ones produced by the Wigner function W = 
(2irÜ)- 1 W{p). 
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Here, v r - s (i) is obtained from the function G (t) substituting the operator Yo 
by r — s in Eq. (|36b|) . 

Since H m ,n is eigenfunction of Y z with eigenvalue i (m + n + 1), the action 
of exp [2Y Z ln G z (i)] on this function produces 

exp [2Y Z \nG z (i)] Iï m ,„ = Gr+" +1 (t) Iï m , n . 

Hence, 

exp ln G z (t)] exp [G (í) Y_] w = ]T ]T A-+J>+i<±* +1 (t) 



, X- 8 ft) / (r + i)!( S + i)! 
X ^^V rïsl "'- 

where u r — s (t) is obtained from G z (í) substituting the operator Yó by r — s in 
Eq. (|36a|) . 

Finally, the action of the exponencial exp [G (i) YL] on the above result gives 



X ÜjWÜl nr+l < s+l ■ 

For the initial condition (|3~4")1 . the Fokker-Planck equation (f2"T)) for the diffusive 
AHO has the following solution 

= T,J2pr+j,s +j K± S s +1 (t)v l r + _Í(t) (37) 
3,1 r > s 



v /( T ·+j)!( s + j)!(r + /)!(7 + 7)! TT 

x Hjw^ Ur+l ' s+l ■ 

Compare this solution with the corresponding one obtained for the density oper- 
ator, Eq. (fTÏÏ)) . or for the Wigner function, Eq. (j2"2")l . in the quantum mechanical 
version of the problem. Mutatis mutandis, the form of the solutions is identical. 
Substituting the functions 7 and Ç in Eq. (|22|) by v and u, we find the result 
given in Eq. (|37|) . 



4 An example 

In order to gain some insight into the differences between quantum and classical 
dynamics of the nonlinear oscillator, let us compare the time evolution of a 
common initial condition. Consider a coherent state |ao)(ao|. The associated 
Wigner function in the complex phase space coincides with the corresponding 
classical distribution, i.e., a gaussian with variance equal to the unity centered 
in «q. An example is shown in Fig. [T] The other figures represent the time 
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Re(a) 



Figure 1: Surface plot of the initial distributions W (a, 0) and w (a, 0). In 
the quantum case, this distribution corresponds to an initial coherent state 
|<2o)(c*o| with ao = 3, whereas in the classical case it corresponds to a gaussian 
distribution centered in a = 3 with unitary variance. 

evolution of this initial state in the quantum and classical models and for the 
regimes with and without difhision. It is important to mention that the time 
evolutions were obtained with the solutions (|22fl and (|37|), 

The classical hamiltonian evolution is such that any point of the phase space 
moves around the origin with angular frequency proportional to |a| , where a is 
the coordinate of this point [5]. Therefore, points over the the initial distribution 
will rotate with an angular velocity that depends on their distance to the origin. 
As consequence, the distribution will continuously spiral around the origin, as 
shown in Fig. [21 The distribution yields a fme-structure in phase space [8], 
which is gradually destroyed if diffusion is included (see Fig. [3]) . 

In the quantum version of the model, the unitary evolution of the Wigner 
function exhibits a very different behavior that the classical one. For times 
equal to rrnr/ (2c/), where m is an integer, the nonlinearity leads to the quantum 
superpositions of states (m is odd) or revivals and anti-revivals (m is even). 
These effects were already reported by Yurke and Stoler [7] and examples of 
them are given in Figs. 2] and [5] 

Whcn the diffusion is included, the quantum effects discussed above are grad- 
ually suppressed. The interference in phase space is reduced, and the regions 
where the Wigner function is negative diminish (see Fig. [6]). For later times, 
the Wigner function and the classical distribution take a form of an annular vol- 
ume around the origin of the phase space. The annular region grows with time 
but its maximum value diminishes in order to mantain constant the integral of 
W (a, t) or w (a, t) over the phase space. This suggests that the Wigner func- 
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-6 -6 

Figure 2: Surface plot of the classical distribution function for the anharmonic 
oscillator at t = 7r(2g) _1 corresponding to the hamiltonian evolution [k = 0) 
of the initial condition shown in Fig. [TJ For this case, g/ui = 0.1. Note the 
fïne-structure yielded by the continuous spiraling of the distribution around the 
origin of the phase space. 

tion of the quantum diffusive AHO converges gradually to the corresponding 
classical distribution function. Non-unitary effects due to quantum dynamics 
of the open AHO were investigated by Milburn and Holmes [9], and by Daniel 
and Milburn [10j . considering the coupling to a null and non-null temperature 
reservoir, respectively. In both works, the authors evaluates the time evolution 
of the Husimi function (Q function), another phase space representation of the 
density operator. Their results show a similar behavior for the Q function, in 
qualitative agreement with the one reported here. 

These results serve to illustrate the procedure and to show that the results 
are physically consistent. However, a detailed study of the quantum to classical 
transition in the AHO must take into account the role played by the parameters 
of interest, namely the nonlinearity strengh, the diffusion constant, and the 
amplitude of the initial coherent state (a measure of classicallity of the initial 
state). This work is in progress. 
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Figure 3: Surface plot of the classical distribution functions for the anharmonic 
oscillator at instants t = n (2g)~ 1 (top) and t = ng' 1 (bottom) corresponding 
to the diffusive evolution (n/g — 0.1) of the initial condition shown in Fig. [TJ 
For this case, g/ui = 0.1. Note that the fine-structure that should be yielded in 
the difïusionless regime (see Fig. [2]) is gradually destroyed. 
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Figure 4: Surface plot of the Wigner function for the quantum anharmonic 
oscillator at t = tt (2g)^ 1 corresponding to the unitary evolution (k = 0) of 
the initial condition shown in Fig. [T] For this case, g/u> = 0.1. Note that 
the appearance of phase space interference due to the coherent superposition of 
states (Schròdinger cat state). At this time, the state exhibits some squeezing 
as well |10j . Compare with the corresponding classical one shown in Fig. [2l 




Figure 5: Surface plot of the Wigner function for the quantum anharmonic 
oscillator at t = ng^ 1 corresponding to the unitary evolution (k — 0) of the 
initial condition shown in Fig. [TJ For this case, g/ui = 0.1. At this instant, 
the oscillator exhibits an anti-revival since it is found in a coherent state with 
amplitude — ao, where ao is the amplitude of the initial coherent state. 
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Figure 6: Surface plot of the Wigner function for the quantum anharmonic 
oscillator at instants t = tt (2g) 1 (top) and t — ng^ 1 (bottom) corresponding 
to the diffusive evolution (n/g = 0.1) of the initial condition shown in Fig. [TJ 
For this case, g/u> = 0.1. Comparing with Figs. H] and [5] one notes that the 
phase space interferences and the revivals are gradually suppressed. 
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